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ABSTRACT 

We identify an abundant population of extreme emission line galaxies (EELGs) at redshift z ^ 
1.7 in the Cosmic Assembly Near-IR Deep Extragalactic Legacy Survey (CANDELS) imaging from 
Hubble Space Telescope/ Wide Field Camera 3 (HST/WFC3). 69 EELG candidates are selected by the 
large contribution of exceptionally bright emission lines to their near-infrared broad-band magnitudes. 
Supported by spectroscopic confirmation of strong [OIII] emission lines - with rest-frame equivalent 
widths ~ lOOOA - in the four candidates that have HST/WFC3 grism observations, we conclude 
that these objects are galaxies with ~ 10* Mq in stellar mass, undergoing an enormous starburst 
phase with M^jM^ of only ~ 15 Myr. These bursts may cause outflows that are strong enough to 
produce cored dark matter proflles in low-mass galaxies. The individual star formation rates and the 
co-moving number density (3.7 x 10^"' Mpc^"^) can produce in ~4 Gyr much of the stellar mass density 
that is presently contained in 10^ — 10^ Mq dwarf galaxies. Therefore, our observations provide a 
strong indication that many or even most of the stars in present-day dwarf galaxies formed in strong, 
short-lived bursts, mostly at z > 1. 
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The formation history of dwarf galaxies with masses '^ 
10* Mq can usually only be studied through 'archaeolog- 
ical' age reconstruction, based on resolved stellar popula- 
tions (e.g., Grebcl 1997; Mateo 1998; Weisz et al. 2011). 
Their high-redshift progenitors have so far remained elu- 
sive despite the ever increasing depth of spectroscopic 
observing campaigns and imaging from the ground and 
the Hubble Space Telescope (HST). In this paper we iden- 
tify an abundant population of z > 1 dwarf galaxies un- 
dergoing extreme starbursts, through HST/ Wide Field 
Camera 3 (WFC3) imaging from the Cosmic Assem- 
bly Ne ar-IR Deep Extragalactic Legacy Survey ( CAN- 
DELS. iGrogin et al.ll20Tlt IKoekemoer et al.|[20lTI ). that 
may well be the progenitors of present-day dwarf galaxies 
with stellar masses ~ 10* — 10^ Mq. 

At the present day, starbursts contribute a minority 
to the total sta r formation activity in dwarf galaxies 
(jLee et al.l 120091 ). However, there is abundant evidence 
that the star formation histories are complex a nd that 
bursts play an important role (as reviewed by IMateol 
1998). Many authors find evidence for short-lived (^^10 
Myr) SF events in nearby star-forming dwarf galaxies 
from a range of observational and modeling techniques 
/e.g.. ISchaereret all 119991: IMas-Hesse fcKuiithI 119991: 
iTho rnlev et al.l 120001: iTremonti et al.l 120011: [Harris et al.l 
.2004 ) , while others argue that star forma tion epochs of 
dwarf galaxies are more prolonged (e.g., iCalzetti et al.l 
ll997trLee 2008 : McOuinn et al. 2009 ). Simulations also 
indicate that star formation histories of low-mass galax- 
ies are episodic or even burst-like (es-- iPelupessv et al.l 
I2OOI IStinson et al.ll2007t INagaminell20m 

As most st ars in dwarf galaxies formed more than 5 
Gyr ago fe.g.. lDolph in et al. 2005; Weisz et al. 2011), it 
is crucial to understand the mode of star formation in 
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Figure 1. Observed I — J vs. J — H colors (in AB magnitudes) 
from HST/WFC3 and ACS imaging for all objects in the UDS 
and GSD (small black points) and the sample of emission-line 
dominated objects (large red points with error bars), selected by 
I - J > 0.44 + (t(7 - J) and J - H < -0.44 - a{J - H), where 
a{I — J) and a{J — H) are the Itr uncertainties on the colors. The 
blue line represents the reds hifted (z = 1.7) continuum colors of 
the Starburst99 model (iLeitherer et al.lll99911 for continuous star 
formation, in the age range from 1 Myr to 100 Myr. The red 
line represents the same model, but with the J-band flux density 
increased by the emission line luminosity predicted by the model 
(Starburst99 predicts Ha luminosity — [OIII] emission, which falls 
in the J band at z = 1.7, is assumed to have the same equivalent 
width). The black arrow indicates dust attenuation. 

dwarf galaxies at those early epochs, but 'archaeologi- 
cal' studies do not have the resolution in terms of stellar 
population age to constrain strengths, durations, and fre- 
quency of bursts. The increased frequency of interactions 
with other galaxies and higher gas fractions at z > 1 may 
have resulted in strong, short-lived starbursts. In this 
paper we place the first constraints on the open ques- 
tion of how many and how frequently strong, short-lived 
starbursts occur in dwarf galaxies at z > 1, and how rel- 
evant this mode of star formation is for the build-up of 
the dwarf galaxy population in a cosmological context. 

2. DATA 

2.1. Multi-Wavelength Imaging 

We select objects from multi- wavelength photometry 
of two fields with HST/WFC3 and Advanced Camera 
for Surveys (ACS) coverage: the Ultra Deep Survey 
(UDS) field and the GOODS-South Deep (GSD) field 
at 4-epoch depthH For the UDS we use WFC3 imag- 
ing in F125W (J) and F160W (H) and ACS imaging in 
F814W (/) from CANDELS. For the GSD we use the J 
and H band imaging from CANDELS, supplemented by 

^^ GOODS is the Great Observatories Origins Deep Sur- 
vey. CANDELS provides deep images over the central parts 
of GOODS-North and GOODS-South, and wider, less deep 
imaging over the remainder of those fields and over the other 
CANDELS fields, including UDS. See the CANDELS website, 
[http://candels.ucolick.org/ , for details of the field layouts. 



WFC3 imaging fro m the Early Release Science (ERS) 
program (Windhor st et al.l 1201 Ih, and / band imaging 
from GOODS (Giavalisco et al.l 120041 ). The total area 
with /-, J-, and H-hand coverage used here is 279 square 
arcminutes. 

Sources are detected i n the H band with SExtractor 
(jBertin fc Arnouti I1996D and p hotometry is performed 
with TFIT (jLaidler et'alll2007[ ). which uses additional 
imaging data sets, ranging from U to 4.5/im to produce 
resolution-matched, multi-wavelength catalogs. The cat- 
alog construction is described in full by Guo et al. (in 
prep.). In addi t ion, w e use a version of GALAPAGOS 
(IHaussler et all 120071) adapted for CANDELS WFC3 
imaging to measure structural parameters (van der Wei 
et al., in prep.). 

2.2. Color-Color Selection 

We select objects that are red in / — J and blue in 
J — H (see Figure [1]), tracing luminous emission lines 
that contribute significantly to the total J-band light. No 
known continuum emission can produce such broad-band 
colors. The highlighted objects in Figure [1] have I — J > 
0.44 -f a{I - J) and J -H < -0.44 - a{J - H), where 
a refers to the color uncertainty; that is, we select those 
objects that are significantly more than 50% brighter in 
J than in both I and H. We identify 69 such objects, that 
is, there is 1 per ~4 square arcminutes. They range in 
magnitude from Hab = 24 to Hab = 27, with a median 
of i? — 25.8 (see Table 1). We note that there is no 
gap in color-color space between the emission line galaxy 
candidates that we select and the general distribution; 
the selected objects are merely the most extreme outliers. 

In Figure [2] we show false-color composites of all 69 
candidates. These sources are typically compact, but 
not unresolved; their J- and 7J-band half-light radii from 
GALFIT are typically 0.1". A subset (- 20%) are more 
extended or consist of multiple components. We show the 
U through A.bfim spectral energy distributions (SEDs) of 
a subset of the emission-line candidates in Figure |3l The 
SEDs are seen to be almost entirely flat in F^, or in terms 
of ultra-violet spectral slope they have /3 ~ — 2, where /3 
is defined as F\ = A^. The J band is a notable outlier 
from this SED shape for all these objects. 

3. EXTREMELY BRIGHT EMISSION LINES 
3.1. Photometric Constraints 

No known objects have continuum SEDs that resemble 
those shown in Figure [31 in particular, the extraordinar- 
ily blue J — H colors are difficult to explain by any radia- 
tive process. Our hypothesis is that the J-band excess 
is due to one or more emission lines. The implied equiv- 
alent widths in the observed frame are extraordinarily 
high: EW ~ 1500 - 3000A. 

Among the emission lines that can reach such extreme 
EWs, Lya and [Oil] are immediately ruled out because 
the implied high redshift would produce a strong break in 
the SEDs; the lack of such a break impli es z < 2 for these 
objec ts. WFC3/UVIS observations ([Windhorst et al.l 
[201 If ) provide UV photometry over the ERS area. The 
average color of those candidates in the ERS area is 
F275W — U ^ 1.44, which suggests that the Lyman 
break is situated at around 3000A in the observed frame, 
which, in in combination with the very blue continuum 
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Figure 2. False-color composites, created from HST /, J, and H band image cutouts of the 69 emission line galaxy candidates. The 
cutouts are 3" on a side, the pixel scale is 0.06", and the full-width half maximum resolution at the longest wavelength (the H band) is 
~ 0.18". The IDs correspond to those in Tabic 1. The sources are typically compact, although a subset of about 20% have more extended 
morphologies or feature multiple components. 



slopes redward of the U band, implies z > 1.5. Only one 
candidate with UVIS coverage has a F275W — U color 
consistent with that of a galaxy at z < 1.5. The implica- 
tion is that strong [OIII] emission at 4959A and 5007A 
provide the most plausible explanation for the J band 
excess light. 

If [OIII] is responsible for the J band excess the redshift 
upper limit is z = 1.8. Furthermore, because we select 
objects with blue J — H colors, the H cannot contain 
the bright Ha line, which implies z > 1.6. Thus, solely 
based on their photometric properties, we suggest that 
our candidates are strong [OIII] emitters in the redshift 
range 1.6 < z < 1.8. 

3.2. Spectroscopic Constraints 

The hypothesis that [OIII] emission at z ^ 1.7 ex- 
plains the J band excess light is strongly supported by 



spectroscopic observations. While none of the candi- 
dates have ground-based spectra, WFC3 grism observa- 
tions are available for sm all portions of the G SD (one 
pointing in the ERS field, IStraughn et al.|[2011[ ) and the 
UDS (from the supernova follow-up program 12099, PI 
A. Ricss). The available grism coverage overlaps with the 
positions of 4 candidates in our sample (1 in the ERS, 3 
in the UDS), and strong emission lines are detected in all 
4 cases. The spectra (Figured]) all show bright emission 
lines in the wavelength range 1.3 — lAjj,m, whose com- 
bined fluxes are in agreement with the excess light seen 
in the J band. 

The lines in all 4 spectra are readily identified as [OIII] : 
the asymmetry of the bright line, always extended blue- 
ward, is due to the two components of the [OIII] line, 
at 5007A and at 4959A, where the latter is 3x fainter. 
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Figure 3. Broad-band SEDs of the 12 emission line galaxy can- 
didates selected from the GSD field. Units on the y-axis are ar- 
bitrary, and the SEDs are incrementally offset by 0.4 dex in the 
vertical direction for clarity, sorted by continuum slope, indicated 
by the dotted lines. The objects are characterized by flat SEDs 
in Fi, over the entire range from U band to H band. The J band 
noticeably deviates from this trend as the result of strong emission 
line contributions, mostly [OIII] at z ^ 1.7. The observed wave- 
lengths for 2 = 1.7 of various emission lines are indicated by the 
vertical dashed lines. 

In all cases, H/3 is also detected. The redshifts are all 
in the range z = 1.65 — 1.80, in excellent agreement 
with what we inferred solely from photometry. We con- 
clude that our sample of extreme emission line galaxies 
(EELGs) form the high-EW tail of the general population 
of emission line galaxies seen in ACS and WFC3 spec- 
troscopic grism observations Ce.g.. [Straughn et al.l 120081 
[200l lAtek et al.ll2010l: IStrauehn et al.ll201lD . 

In principle, Ha emitters at 0.9 < z < 1.1 should also 
be included by our selection technique. The spectroscopy 
and UV photometry indicate that those must be far fewer 
than [OIII] emitters. Whether this is due to selection ef- 
fects or evolution in the number density of such objects 
remains to be seen. Still, even though our working hy- 
pothesis is that all 69 candidates are [OIII] emitters at 
z ~ 1.7, we should keep in mind that some fraction of 
our 69 candidates are likely Ha emitters at z ~ 1. 

3.3. Emission Lines and Broad-Band Photometry 

We have shown that selecting objects which are much 
brighter in J than in / and H works as a rather clean 
method for finding strong [OIII] emitters at 1.6 < z < 
1.8. Emission line galaxies with such excesses in other 
bands also exist, but a systematic search is more compli- 
cated as at most redshift ranges, multiple lines (most no- 
tably [OIII] and Ha) affect multiple photometric bands. 
Therefore, we refrain from conducting such a systematic 
search here. 

The existence of such emission-line dominated galax- 
ies complicates the interpretation of SEDs, which is es- 
pecially relevant in the case of the search for and SED 
modeling of rare, high-redshift objects. Although con- 
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Figure 4. WFC3 grism spectra of the four candidates with grism 
coverage. The IDs refer to those in Table 1. GSD18 is object 402 
from Strauglm et al. (2011); the 3 objects in the UDS are from 
supernova follow-up grism observations (program ID 12099, PI 
A. Riess). The three vertical dashed lines show positions of H/3, 
[OIII], and Ha for z = 1.7. These spectra strongly suggest that the 
majority of the objects in our sample are [OIII] emitters at 2 ~ 1.7. 

tamination by emission lines is often considered to be 
a factor (e.g., Labbc et al. 2010), the extremely bright 
lines we obser ye suggests that their effect may be un- 
derestimated. lOno et all ()2010[ ) explicitly showed that 
red colors in Lya emitters and z — 7 Lyman break 
galaxies may indicate the presence of evolved stellar 
populations or strong nebular emission lin es (also see 
iSchaerer et"allll999l : iFinkelstein et al.ll2011f ). Steep UV 
continuum slopes, such as observed in our objects, should 
serve as a warning sign for contamination by nebular 
emission lines at longer wavelengths to the point that 
those can dominate the broad-band fiux density. 

4. STARBURSTING DWARF GALAXIES AT Z = 1.7 
4.1. Star formation or AGN? 

Before turning to our preferred starburst interpreta- 
tion, let us first point out that nuclear activity is not a 
likely explanation for the bright emission lines in the vast 
majority of EELGs. None of the objects in the CDFS 
have significant detections in X-ray or at 24/im. The ob- 
jects are spatially resolved in both J and H, and, more- 
over, the J and H band sizes are consistent with each 
other. 

Moreover, it is highly unlikely that all 69 objects 
are dominated by line emission from active galactic 
nuclei. At least at the present day, low-mass, low- 
metallicity AGN are exceedingly rare (jlzotov fc ThuanI 
'2008^, much rare r than starbursting dwarf galaxies 
(Izo tov et al.ll2011l ). The implied black hole masses for 
the objects in our sar nple, as inferred f rom their UV con- 
tinuum luminosities ()Shen et al.ll2008t ) are ^ 10^ Mq at 
most, when assuming an Eddington accretion of unity. 
At these low masses, at least at the present day, secular 
processes drive nuclear activity; thus, an unknown accre- 
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Figure 5. Rest-frame V-band absolute magnitude vs. rest-frame 
equivalent width of the [OIII] emission line at 5007A as calculated 
from the broad-band photometry as explained in the text and as- 
suming that all emission-line galaxy candidates are at 2 = 1.7. 
They span a range in luminosity, My = —17 to My = —20, and 
have EW[o 1 1 1],5007 between 500A and 1200A. See also Table 1. 

tion mode or triggering mechanism for nuclear activity 
would have to be invoked in order to explain an extreme 
change in the relative numbers of AGN and starburst 
powered emission-line dominated objects. At these low 
masses, merging cannot account for this. The starburst 
hypothesis, on the other hand, places these objects in the 
realm of dwarf galaxies, and their physical and statisti- 
cal properties are consistent with the abundances and 
masses of dwarf galaxies as we will discuss below. 

Although nuclear activity cannot be ruled out entirely 
- and line-strength gradients in star forming z ^ 2 galax- 
ies s uggest that weak A GN may contribute to some ex- 
tent ([Trump et al.ll20fll ) - we can safely assume that the 
observed emission lines are effectively dominated by star 
formation activity. 

4.2. Starburst Ages and Masses 

We interpret the observatio ns in the contex t of the 
Starburst99 model (SB99, Leit herer et"aI]|1999D . which 
includes predictions for how the EWs of Hydrogen re- 
combination line evolve over time. Therefore, our first 
task is to estimate H/3 line strengths from the data. We 
attribute the excess light in the J band, compared to the 
continuum light measured in the I and H bands, to com- 
bined effect of emission lines. Therefore, we can compute 
the combined equivalent width as follows: 



EW 



f.T 



fr + fF 



Wj 

1 + z 



(1) 



where Wj = 2845A is the effective width of the J-filter 
response curve, z = 1.7 to correct the observed EW to the 
rest frame, and / is the flux density f„ in the respective 
filters. 
The relative contributions of the various emission lines 



are constrained by fitting Gaussian components to the 3 
emission lines seen in the grism spectra shown in Figure 
|4j keeping the ratio between the two [OIII] components 
fixed at 3. We only use the 3 UDS spectra as H/3 is only 
marginally detected in the GSD spectrum. The emission 
line ratios are remarkably similar for all 3 objects: II/3 
contributes 1/8 to the combined line l uminosity, suggest- 
ing a very low met allicity (see, e.g., iSalzer et al.l 120051 : 
lAmorin et al.l [20Tol for comparisons). Because the fiux 
is dominated by the [OIII] 5007 line and is therefore more 
directly related to our observations, we show the inferred 
[OIII] 5007 EWs in Figure [5] (also see Table [J). How- 
ever, we model the observations by fitting the inferred 
H/3 EWs to the SB99 predictions. These are assumed to 
be always l/8th of the combined EW. The unavoidable 
intrinsic scatter in this conversion is mimicked by prop- 
agating a generous factor of two in the uncertainties of 
the quantities we infer below. 

EW-ap is a sensitive age indicator, as it is quickly re- 
duced once a stellar population gains in mass or the star 
formation activity diminishes. F or a SB99 model with 
continuous star formation with a iChabrien (|2003| ) IMF 
with a high-mass cut off at IOOM0 and metallicity 0.2Zq 
the H/3 EWs imply that the galaxies in our sample typ- 
ically have ages of 10 — 20 Myr (Figure [S]). It seems 
unlikely that these bursts will be much longer than this 
given the energy put into the interstellar medium (see 
Sec. 4.3); follow-up grism observations will directly con- 
strain the number of older bursts. 

If we assume a single burst model instead, we infer 
ages 3 — 5 Myr; all formation histories with declining 
star formation rates produce ages that are bracketed by 
these two extremes. Continuous star formation seems 
more realistic than an instantaneous burst, a notion that 
becomes physically untenable at very young ages. Star 
formation cannot happen faster than the dynamical time 
scale of '^lO Myr for these systems. Most relevant for our 
analysis is that the stellar mass estimates derived from 
the two different models are very similar (see below). 

We note that choosing a different metallicity does not 
significantly change our results - a low metallicity is re - 
alistic for these low-mass systems (|Amorin et al.l 120101 ). 
and the [0III]-H/3 ratio suggests the metallicity is indeed 
low. Choosing an even lower metallicity (0.05Zq instead 
of 0.2Zq) results in an increase in both age and mass by 
'^0.2 dex. We do suffer from the usual, unknown uncer- 
tainty due to our lack of knowledge of the stellar IMF. A 
different slope or cutoff at the high-mass end changes the 
number of ionizing photons from high- mass stars per unit 
stellar mass, and the entirely unconstrained number of 
low-mass stars (< 1 Mq) determines the overall normal- 
ization of the stellar mass. In general, we conclude that 
the ages of these galaxies are < 40 Myr, which includes 
the intrinsic range in age and the systematic uncertainty 
due to the unknown star formation history. 

In the following we use the results from the continuous 
star formation model, but using the instantaneous burst 
model, by virtue of the insensitivity of the mass estimates 
to the choice of star formation history, does not change 
our interpretation and conclusions. 

For a given age, the SB99 model predicts the (rest- 
frame) y-band mass-to-light ratio, such that we can di- 
rectly estimate the mass after deriving the T^-band lumi- 
nosity from the observed H band magnitude. We cor- 
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Figure 6. Masses, ages and star formation rates for the 69 
emission-line dominated objects in our sample, derived with the 
SB99 model, assuming that all are at z = 1.7, adopting a con- 
tinuous SF model with 0.2 times solar metallicity and a Chabrier 
IMF. The SFRs indicated by the diagonal lines are simply obtained 
by dividing the mass (x-axis) by the age (y-axis). The galaxies 
in our sample typically have 10* Mq stellar masses with young 
ages (5 — 30 Myr), or, equivalently, extremely high specific SFRs 



(~ 5 X IQ-* yr-l, or ~ 50 X t\ 



J- 



rect the luminosity and the derived mass estimate for 
extinction by comparing the continuum slope derived 
from the ACS photometry at rest-frame 2500A, typi- 
cally /32500 2 (see Table 1), with the SB99 model 

prediction (rather constant at &if, m "^ — 2.6 for the ages 
of these bursts). If we adopt the iCalzetti et all ()2000l ) 
extinction law for starbursting galaxies the typical ex- 
tinction is E{B -V)^ 0.2. 

The median mass we infer is 8 x 10^ M0 (see Table 
1 and Figure H]) . Mass estimates inferred from the in- 
stantaneous burst model are only slightly smaller, by 
less than 0.1 dex on average. Internal consistency lends 
our modeling approach strong credibility: given the in- 
ferred ages, masses, and extinction corrections, the SB99 
model predicts dust- attenuated rest-frame UV luminosi- 
ties that are consistent with the observed rest-frame UV 
luminosities - the latter are not used in our modeling 
procedure. Thus, the model successfully describes the 
observed rest-frame UV and optical continuum spectral 
energy distributions as well as the observed emission line 
luminosities. 

Full modeling of the spectral energy distributions that 
includes emission line contributions will be presented in 
forthcoming studies that will also include objects with 
less prominent emission lines. As a consistency check 
with the results presen ted above we al r eady applied 
the method outlined by iFinkelstein et al.1 (|2011[ ) to the 
galaxies in our sample. This method compares the ob- 
served photometry with model spectral energy distribu- 
tions that include contributions from nebular and stel- 
lar continuum radiation as well as all nebular emission 



lines. Free parameters in this modeling procedure in- 
clude redshift, stellar mass, extinction and metallicity. 
We find photometric redshifts that are consistent with 
1.6 < z < 1.8 for the vast majority of the sample. More- 
over, the inferred stellar masses and ages are very simi- 
lar, even though this model is based on a different stellar 
population synthesis model. 

The implied star-formation rates (see Figure |6]) may 
lead one to expect that these galaxies have significant 
24/xm detections. However, we have verified that none of 
the 29 candidates in the GSD, which has very deep 24/im 
from GOODS, is detected. A possible explanation for 
this is the presumably low metallicity, which would result 
is relatively small dust masses and, hence, low infrared 
luminosities. 

The observed IRAC flux densities at 3.6/ini and 4.5/im 
from GOODS and SED^3 are in most cases - there are 
two exceptions - fully consistent with the expected flux 
densities for the bursts observed in the UV. In addition, 
the galaxies have the same sizes in the J and H bands, in- 
dicating that the spatial extent of the region from which 
the line emission originates roughly follows the stellar 
light. Hence, there is no evidence for underlying older 
stellar populations. However, we cannot rule out their 
existence: maximally old stellar populations have mass- 
to-light ratios that are up to ~ 50 times larger than those 
of the bursts, even in the near infrared. 

If we assume a past star formation rate that is con- 
stant after averaging over > 100 Myr time scales we find 
upper limits for the mass in older stars that is ~ 5x the 
burst mass. The implied total stellar mass upper limits 
are then ^ 5 x 10^ Mq. This caveat notwithstanding, 
we assume in the remainder of this paper that there is 
no significant population of older stars in these galaxies, 
and that the observed bursts account for the total stellar 
mass. However, the bottom line is that the total stellar 
masses of these objects are well below 10^ Mq, in the 
regime of dwarf galaxies. 

5. DISCUSSION 

5.1. Comparison with Other Samples 

Galaxies with similar properties have previously been 
identified through broad-band p hotometry at z < 0.4 i n 
the Sloan Digital Sky Survey (Cardamone et al. '2009), 
and have been s hown by Amori'n et al. (2010) and 
llzotov et al.l (j2011f) to constitute the most strongly star- 
forming and most metal poor tail of th e well-known class 
of blue compact dwarf galaxies fe.g.. iSargent fc Searld 
[T970I: iTblan fc Martinl[l98ll: IGriffith et al.ll201lD . which 
have very low metallicities and extremely high, spatially 
concentrat ed star-format ion activity (jGuzman et al.l 
11998: Overzier et al.ll2008l) . 

iCowie et al.l (|2011f ) (also see iScarlata et aH [2009)) 
recently studied the Lya properties of high-EW Ha 
emitters, providing a direct connection between higher- 
redshift searches of Lya (e.g.. lOuchi et al.|[2008l : lHu et al.l 
I2010( ). and find Lya EWs rangin g from 20A to 200A . 
Combining this with the findings of lNilsson et al.l (|2011| ) , 
who show that Lya emitters at z ~ 2 are objects with a 
very wide range in properties, it is clear that from Lya 

21 GOODS and SEDS (Spitzer Extended Deep Survey) provide 
the deepest IRAC imaging ever obtained, an excellent probe of 
stellar populations at high redshift 
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Table 1 

Sample of Extreme Emission Line Galaxies 



EELG2011 



RA (J2000) 
(deg) 



DEC (J2000) 
(deg) 



H 

(AB) 



-£'W'[OIII],5007 

(A) 



/32500 



log (Af ) 



GSDl 


53.167064 


-27.858936 


24.67 ± 0.07 


459 ± 40 


-1.83 ± 0.09 


8.57 ± 0.17 


GSD2 


53.080345 


-27.850572 


25.63 ± 0.07 


569 ± 67 


-2.04 ± 0.19 


7.96 ± 0.19 


GSD3 


53.046020 


-27.837322 


25.88 ± 0.11 


507 ± 75 


-2.02 ± 0.18 


7.94 ± 0.20 


GSD4 


53.105087 


-27.819974 


26.20 ± 0.11 


769 ± 143 


-1.75 ± 0.30 


7.72 ± 0.22 


GSD5 


53.067508 


-27.773595 


25.05 ± 0.08 


566 ± 74 


-2.19 ± 0.19 


8.14 ± 0.19 


GSD6 


53.097499 


-27.763919 


24.99 ± 0.05 


700 ± 53 


-2.12 ± 0.09 


8.10 ± 0.17 


GSD7 


53.122127 


-27.759542 


25.44 ± 0.12 


535 ± 99 


-1.70 ± 0.28 


8.20 ± 0.22 


GSD8 


53.171936 


-27.759145 


24.26 ± 0.04 


693 ± 47 


-1.76 ± 0.11 


8.52 ± 0.17 


GSD9 


53.078754 


-27.750288 


24.86 ± 0.04 


468 ± 32 


-1.99 ± 0.10 


8.42 ± 0.17 


GSDIO 


53.063690 


-27.745853 


26.37 ± 0.09 


759 ± 134 


-1.56 ± 0.33 


7.72 ± 0.22 


GSDll 


53.007499 


-27.741867 


25.97 ± 0.09 


534 ± 76 


-2.13 ± 0.19 


7.84 ± 0.20 


GSD12 


53.114612 


-27.721979 


25.80 ± 0.12 


641 ± 139 


-1.94 ± 0.31 


7.87 ± 0.23 


GSD13 


53.101516 


-27.720882 


24.77 ± 0.03 


490 ± 29 


-2.35 ± 0.08 


8.29 ± 0.16 


GSD14 


53.055908 


-27.718803 


26.00 ± 0.08 


501 ± 65 


-1.90 ± 0.20 


7.94 ± 0.19 


GSD15 


53.149536 


-27.710285 


26.64 ± 0.15 


820 ± 288 


-2.20 ± 0.57 


7.36 ± 0.30 


GSD16 


53.147617 


-27.707088 


26.10 ± 0.08 


582 ± 80 


-2.08 ± 0.22 


7.75 ± 0.20 


GSD17 


53.064220 


-27.706523 


25.41 ± 0.05 


465 ± 43 


-2.20 ± 0.13 


8.12 ± 0.17 


GSD18 


53.071292 


-27.705802 


25.24 ± 0.04 


861 ± 66 


-2.36 ± 0.11 


7.85 ± 0.17 


GSD19 


53.181976 


-27.705038 


25.71 ± 0.06 


1002 ± 245 


-2.18 ± 0.41 


7.72 ± 0.23 


GSD20 


53.140815 


-27.692390 


26.23 ± 0.10 


496 ± 82 


-1.83 ± 0.26 


7.88 ± 0.21 


GSD21 


53.100936 


-27.676704 


24.76 ± 0.10 


935 ± 139 


-0.95 ± 0.28 


8.54 ± 0.21 


GSD22 


53.118450 


-27.819919 


26.76 ± 0.13 


870 ± 198 


-1.88 ± 0.36 


7.41 ± 0.24 


GSD23 


53.077606 


-27.812795 


26.81 ± 0.17 


1512 ± 338 


-2.22 ± 0.28 


7.30 ± 0.25 


GSD24 


53.132972 


-27.740102 


27.77 ± 0.30 


698 ± 318 


-2.23 ± 0.47 


6.95 ± 0.38 


GSD25 


53.084388 


-27.727920 


27.29 ± 0.15 


562 ± 164 


-2.39 ± 0.43 


7.18 ± 0.27 


GSD26 


53.141502 


-27.724880 


26.65 ± 0.10 


650 ± 106 


-2.10 ± 0.25 


7.47 ± 0.21 


GSD27 


53.112579 


-27.707090 


26.91 ± 0.13 


954 ± 262 


-2.40 ± 0.42 


7.15 ± 0.26 


GSD28 


53.046119 


-27.705604 


27.13 ± 0.16 


1009 ± 293 


-2.14 ± 0.45 


7.16 ± 0.28 


GSD29 


53.139953 


-27.675138 


27.79 ± 0.21 


1314 ± 557 


-2.12 ± 0.74 


6.93 ± 0.37 


UDSl 


34.275299 


-5.274496 


25.38 ± 0.09 


576 ± 90 


-1.33 ± 0.18 


8.31 ± 0.20 


UDS2 


34.440769 


-5.262566 


25.74 ± 0.09 


1081 ± 147 


-1.41 ± 0.17 


7.99 ± 0.20 


UDS3 


34.482173 


-5.261399 


25.28 ± 0.08 


507 ± 95 


-1.88 ± 0.24 


8.23 ± 0.20 


UDS4 


34.268657 


-5.260064 


25.44 ± 0.10 


614 ± 83 


-1.53 ± 0.14 


8.18 ± 0.19 


UDS5 


34.426483 


-5.255770 


25.69 ± 0.11 


701 ± 95 


-1.66 ± 0.10 


7.99 ± 0.20 


UDS6 


34.428569 


-5.255318 


25.10 ± 0.07 


731 ± 86 


-2.12 ± 0.13 


8.04 ± 0.18 


UDS7 


34.325676 


-5.251743 


24.32 ± 0.04 


656 ± 43 


-1.59 ± 0.05 


8.58 ± 0.17 


UDS8 


34.314014 


-5.251047 


26.44 ± 0.17 


728 ± 153 


-1.39 ± 0.20 


7.77 ± 0.24 


UDS9 


34.382587 


-5.244620 


25.94 ± 0.09 


478 ± 64 


-1.82 ± 0.14 


8.03 ± 0.19 


UDSIO 


34.263534 


-5.239433 


25.48 ± 0.07 


541 ± 64 


-1.84 ± 0.14 


8.13 ± 0.18 


UDSll 


34.311279 


-5.238957 


26.36 ± 0.10 


735 ± 94 


-2.43 ± 0.12 


7.42 ± 0.19 


UDS12 


34.473888 


-5.234232 


24.15 ± 0.03 


713 ± 42 


-1.72 ± 0.05 


8.57 ± 0.16 


UDS13 


34.318141 


-5.232299 


25.35 ± 0.07 


716 ± 68 


-2.12 ± 0.08 


7.95 ± 0.18 


UDS14 


34.481567 


-5.222499 


25.69 ± 0.11 


602 ± 96 


-2.30 ± 0.15 


7.81 ± 0.20 


UDS15 


34.371166 


-5.214803 


25.45 ± 0.09 


843 ± 111 


-1.35 ± 0.15 


8.14 ± 0.19 


UDS16 


34.482921 


-5.214187 


25.39 ± 0.08 


662 ± 87 


-1.80 ± 0.15 


8.07 ± 0.19 


UDS17 


34.247516 


-5.205330 


25.95 ± 0.10 


469 ± 63 


-2.22 ± 0.13 


7.89 ± 0.19 


UDS18 


34.315448 


-5.200902 


25.83 ± 0.13 


739 ± 136 


-1.87 ± 0.17 


7.84 ± 0.21 


UDS19 


34.298866 


-5.191800 


26.25 ± 0.15 


543 ± 106 


-1.94 ± 0.17 


7.79 ± 0.23 


UDS20 


34.232082 


-5.190388 


25.16 ± 0.06 


648 ± 55 


-2.20 ± 0.06 


8.03 ± 0.17 


UDS21 


34.308940 


-5.190090 


25.15 ± 0.07 


609 ± 66 


-2.28 ± 0.11 


8.03 ± 0.18 


UDS22 


34.416740 


-5.180443 


26.34 ± 0.19 


1070 ± 307 


-1.72 ± 0.31 


7.64 ± 0.28 


UDS23 


34.387023 


-5.177240 


25.59 ± 0.07 


591 ± 74 


-2.15 ± 0.14 


7.91 ± 0.18 


UDS24 


34.252845 


-5.176362 


26.02 ± 0.13 


779 ± 140 


-1.40 ± 0.19 


7.91 ± 0.22 


UDS25 


34.402145 


-5.175352 


24.47 ± 0.05 


507 ± 39 


-2.11 ± 0.06 


8.48 ± 0.17 


UDS26 


34.459190 


-5.174448 


25.34 ± 0.08 


552 ± 59 


-2.29 ± 0.09 


8.01 ± 0.18 


UDS27 


34.284194 


-5.164084 


26.40 ± 0.17 


576 ± 139 


-1.05 ± 0.30 


8.01 ± 0.26 


UDS28 


34.500236 


-5.155595 


25.89 ± 0.14 


519 ± 100 


-1.77 ± 0.21 


8.01 ± 0.22 


UDS29 


34.263267 


-5.152174 


26.35 ± 0.13 


1003 ± 150 


-1.90 ± 0.11 


7.57 ± 0.21 


UDS30 


34.477771 


-5.147521 


25.66 ± 0.14 


533 ± 103 


-0.65 ± 0.21 


8.50 ± 0.22 


UDS31 


34.296325 


-5.144416 


25.95 ± 0.12 


546 ± 88 


-2.17 ± 0.14 


7.81 ± 0.20 


UDS32 


34.419242 


-5.142892 


25.83 ± 0.12 


721 ± 138 


-1.69 ± 0.21 


7.91 ± 0.22 


UDS33 


34.246810 


-5.139120 


26.06 ± 0.11 


553 ± 85 


-1.54 ± 0.19 


7.99 ± 0.20 


UDS34 


34.371933 


-5.137272 


24.83 ± 0.05 


586 ± 50 


-2.73 ± 0.08 


8.01 ± 0.17 


UDS35 


34.314693 


-5.133675 


26.42 ± 0.21 


1125 ± 315 


-1.20 ± 0.28 


7.80 ± 0.28 


UDS36 


34.261756 


-5.134672 


25.32 ± 0.08 


658 ± 76 


-2.16 ± 0.12 


7.97 ± 0.18 


UDS37 


34.380569 


-5.268105 


26.54 ± 0.22 


832 ± 249 


-1.34 ± 0.28 


7.71 ± 0.29 


UDS38 


34.441444 


-5.215963 


26.60 ± 0.19 


912 ± 208 


-2.14 ± 0.14 


7.38 ± 0.25 


UDS39 


34.334854 


-5.177163 


26.99 ± 0.18 


594 ± 145 


-2.39 ± 0.24 


7.26 ± 0.26 


UDS40 


34.438095 


-5.160070 


26.86 ± 0.17 


677 ± 157 


-2.29 ± 0.21 


7.30 ± 0.25 



Note. — EELG2011: identification number prefixed by the respective field acronyms; RA/DEC: coordinates 
from the CANDELS catalogs; H: H-haxid AB magnitude from the CANDELS catalog; £'VFroiii],5007' rest- 
frame equivalent width inferred from the /, J, and H broad-band photometry (see text for details); /32500- ^A 
continuum slope at rest-frame 2500A inferred f rom a linear fit to th e S, V, and / broad-band photometry; 
log (A/): stellar mass inferred from Starburst99 dLeitherer et aUI 19991) , as described in the text. 



Extreme Emission Line Galaxies 



emitters one cannot derive a complete description of star 
formation in low-mass galaxies. On the other hand, Lya 
emitters at higher redshifts (z > 3) appear to be young , 
with small stellar masses (e.g., iFinkelstein et al.1 120091 ) . 
similar to the emission line galaxies studied here. 

Narrow-band surveys identified galaxies with strong 
[OIII] and Ha emission line s with EW ^ 100 - lOOOA at 
redshifts z = 0.3 — 1 (e.g., iKakazu et al.ll2007f). demon - 
strated to be young and r netal poor (JHu et al.l [20091 ). 
Most notably, lAtek et all ([2010') pointed out the exis- 
tence of a class of emission line galaxies at z ~ 1.5 with 
EW > lOOOA that would most likely be included in our 
sample as well. However, so far, their nature has not been 
described and their cosmological relevance in the context 
of galaxy formation has remained unclear. Therefore, let 
us now put these starbursting dwarf galaxies in a cosmo- 
logical context. 

5.2. Cosmological Context: Implications for the 
Formation of Dwarf Calaxies 

Our sample with redshifts 1.6 < z; < 1.8 consists of 69 
low-mass (~ 10® Mq), young (^ 0.5—4x10^ yr), extreme 
starbursting, presumably metal-poor galaxies. Their co- 
moving number densit-sFl is 3.7 x lO^'* Mpc~'^, two or- 
ders of magnitude higher than that of n earby galaxies 
with similar EWs (jCardamone et al1l2009| ) . The individ- 
ual star formation rates and the number density com- 
bine into 1.7 X 10-3 Mq yr"! Mpc'^ This is a 2% 
contribution to the total star-formation rate density at 
z ^ 1.7 occurring in galaxies that contribute perhaps 
~ 0. 1% to the total ste llar mass density at that epoch 
(e.g-. lKarimet al.ll20li[) . 

Placing these detections in the context of the burst 
fraction among equally massive dwarf galaxies is diffi- 
cult, however, since we cannot currently constrain their 
number density; galaxies without starbursts have older 
luminosity weighted ages and, as a consequence, are at 
least 2 magnitudes fainter than the starbursting galaxies. 
With marginal detections, even in the latest WFC3 data, 
and no obvious spectral features, their redshifts cannot 
easily be estimated. Model predictions differ strongly 
from observati onal measuremen ts of the galaxy stellar 
mass function (|Guo et al.1 1201 ID and cannot be used to 
constrain the burst fraction among the population of low- 
mass galaxies. 

Nevertheless, we can still gauge the importance of the 
observed burst for the formation of low- mass galaxies by 
making the reasonable (and testable) assumption that 
the observed bursts occur with the same frequency at all 
epochs 1 ^ z < 3.5 - a period of ~4 Gyr during which 
the cosmic star formation history peaked and after which 
the number density of starbursts declines as mentioned 
above. The basic indication that such bursts are im- 
portant is that the number of stars produced in such 
bursts over a period of several Gyr is comparable to the 
number stars in present-day dwarf galaxies. Let us use 
this consideration to construct a simple toy model that 
relates the observations presented here to the mass func - 
tion of present-day low-mass gal a xies. [Guo et al.l (|2011f ) 
use the data from iBaldrv et al.l (|2008f ) for the galaxies 

^^ the values for our two widely separate fields, UDS and GSD, 
differ by only 12% 



with masses down to 10^ Mq, and their mass function 
can be represented by a simple power law for all galaxies 
wit h masses < 10^*^ Mq, that is, well below the knee of 
the lSchechtCTJ (|1976D function: 



'^^''*^-'HJi^) 



-0.37 



(2) 



in units of Mpc"^ dlog(M*)-^ 

Let us then express the stellar mass of the present- 
day descendants of the observed starbursting galaxies at 
z ^ 1.7 in terms of the following star formation history: 



Mdc 



Mburst X Nhu 



/l 



burst 



(3) 



where Afburst --2x15 Myr x 5 M© yr^^ = 1.5 x 10* Mq is 
the total stellar mass produced in a single starburst (the 
factor 2 is included to convert the observed burst age to 
its total duration), iVburst is the number of bursts that 
occurs in each galaxy over the ~4 Gyr period between 
z = 1 and z — 3.5, and /burst is the fraction of the 
total stellar mass that is produced in such bursts - the 
rest is assumed to form in smaller bursts and/or a more 
quiescent mode of star formation. 

Now we can predict the number of present-day de- 
scendants by dividing the number of observed bursts at 
z ~ 1.7 by the duty cycle, that is, the fraction of the 
time that a bursts is occurring over the period of time 
that they can occur (here assumed to be 1 < z < 3.5, or 
4 Gyr). This duty cycle is the ratio of the duration of a 
single burst (2 x 15 Myr) and this 4 Gyr period, multi- 
plied by the number of bursts per galaxy iVburst- Thus 
we have 



4>{Mdcsc) 



0burst 4 Gyr 
iVburst "" 30 Myr' 



(4) 



where (/'burst is the co-moving number den- 
sity of the bursts observed at z ^ 1.7 
(3.7 X 10"^ Mpc~3 dlog(M*)-\ where we use the 
observed span in stellar masses of an order of magnitude 
- see Figure 6 - to introduce this differential unit). 

According to Eqs. 2 and 4 the mass of the descendant, 
Mdcsc, is uniquely determined by the fraction of stars 
formed in bursts, /burst- Equating (p{Mdesc) (Eq. 4) and 
0(M*) (Eq. 2), and substituting for A'^burst (Eq- 3) we 
obtain, after rearranging terms: 



108 M, 



2-4/, 







burst' 



(5) 



Thus, if /burst is close to unity, that is, almost all stars 
are formed in bursts, then we infer that each galaxy must 
undergo one or two bursts on average and that, therefore, 
-^^dosc ~ 1 — 2Mburst- It is perhaps more realistic to adopt 
a smaller value for /burst- If /burst ~ 0.5, we find that 
two or three bursts must occur in each galaxy, producing 
descendants with masses ^ 10^ Mq. The latter would 
imply a growth in stellar mass between z ~ 1.7 and the 
present by at least a factor 3 given the mass constraints 
on the underlying populations of the observed starburst 
galaxies. Galaxy formation models suggest that the typ- 
ical growth is indeed a factor of 3 or 4, but we should 
bear in mind that these models do not reproduce the ob- 
served co-moving number density evolution of low-mass 
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galaxies with redshift (e.g.. lGuo et al.|[2011[ ). Therefore, 
these predictions should be treated with care. 

Choosing /burst very low (< 0.1) implies a very 
large growth in mass, with high-mass descendants (> 
10^° Mq). Such growth by, say, more than an order 
of magnitude in mass seems unlikely in the context of 
current galaxy formation models. In particular, models 
are better observationally constrained for these higher 
masses and the prediction is that more massive galaxies 
j^^Qio-ii Mq) grow in mass by a factor 3 or 4, not by two 
order s of magnitu de, between z ~ 1.7 and the present 
(e.g.. iKeres et al.l [20091. Given these constraints, the 
general conclusion we can draw is that our observations 
suggest that many or most stars in present-day dwarf 
galaxies (with masses < 10^ -^o) have formed in a small 
number of starbursts at z > 1. 

The quantitative interpretation we offer here relies on 
our estimate of the characteristic burst duration (30 
Myr). However, the uncertainty in this assumed burst 
duration does not affect our conclusion that intense star- 
bursts at high redshift constitute an important phase in 
the mass buildup low-mass present-day galaxies. 

If the bursts last significantly longer, which we deem 
unlikely, then the total mass formed in a burst will be 
larger than assumed above, while the frequency of the 
bursts will remain the same since our selection method 
is only sensitive to young bursts. As a result, the fraction 
of stars formed in bursts, /burst will be larger than what 
we derived above. 

Significantly shorter bursts (say, ~ 5 — 10 Myr, as al- 
lowed by the SB99 model - see Sec. 4.2) are also un- 
likely, as the bursts are spatially extended, and the im- 
plied crossing time in these systems is about 30 Myr. 
Nonetheless, in the case that burst times are short, the 
burst frequency must be higher to account for their ob- 
served number (Eq. 4), but the total mass produced in 
each burst remains the same (see Sec. 4.2). As a result, 
the number of bursts per galaxy (A'burst) and, therefore, 
the fraction of stars formed in bursts (/burst) increase. 

The main assumptions in our interpretation, then, are 
1) that the descendants of the observed galaxies remain 
low-mass galaxies up to the present day (< 10^ M© in 
stars), and 2) that the observed bursts do not only oc- 
cur atz = 1.6 — 1.8, but are equally frequent over a 
much broader redshift range (1 < z < 3.5). The first as- 
sumption is supported by the understanding that more 
massive galaxies are not expected to grow in stellar mass 
by several orders of magnitude, as would be required if 
the descendants of the observed starbursting galaxies are 
much more massive. The second assumption is straight- 
forward to test observationally by searches for similar 
objects over a wider redshift range. 

5.3. Energy Budget and Core Formation 

The most remarkable property of these galaxies are 
their growth rates (specific star formation rates) of 20 — 
200 Gyr"""^. This is far outside the realm of normal, more 
massive star forming galaxies, which typically have ~ 
1 Gyr-i ('e.g.. lKarimeran i201ll). 

The amount of energy deposited into the interstellar 
medium through winds and supernovae (10^^~^^ erg) ex- 
ceeds the binding energy by an order of magnitude, even 
if the total mass is larger than the stellar mass by many 



factors. This implies that the gas that is fueling the 
starburst may be in the process of being blown out, re- 
sulting in the end of the starburst phase. Whether the 
gas will be blown out of the halo depends on the halo 
mass, which is currently unconstrained. If the gas does 
not leave the halo it may eventually once more cool and 
sink to the center of the potential well and trigger an- 
other starburst. Then, the observed bursts could be 
part of a semi-periodic cycle of star-formation activity, 
which is seen in simulations of low-mass galaxies (e.g., 
Stinson et al. 2007). Episodic star formation would also 
alleviate the difficulty of explaining the sudden occurence 
of such a large starburst in such small galaxies. If such 
starbursts are not cyclical, then another possible expla- 
nation, usually offered to explain the lack of star forma- 
tion at early times in even less massive systems, is that 
star formation had been suppressed at earlier epo chs as 
a result of UV background rad iation (e.g.. Babul fc ReesI 
ri99l iBabui fc Fer gusoni 119961) . 

An intriguing possibility is that if most of the gas is 
expelled from the central region, then the stars them- 
selves could become unbound as well, dissolving the en- 
tire galaxy. Displacement of large amounts of mate- 
rial leads to changes in the potential, and if this hap- 
pens at very short time scales, then even the dark mat- 
ter profile can be altered. Outfiows have long been ar- 
gued to play a role in prod ucing cored density profiles 
in low-mass galaxies (e.g., iNavarro et al.l 119961 ). The 
time scales, masses, and star-formation rates we derived 
above for the galaxies in our sample roughly meet the 
requirements for such a proc ess to be efficient. Recently, 
iPontzen fc Governatol ()2011|) showed that episodic star 
formation, even with less intense bursts, can produce the 
same effect. What the previous and subsequent evolution 
of the systems we observed is remains to be seen, but the 
energy balance and the possibility that such bursts are 
reoccurring make our observations consistent with this 
general picture; indeed outflows may be responsible for 
the formation of flattened dark matter proflles in low 
mass galaxies. These speculative scenarios can be tested 
further with observational constraints on the gas masses 
and hydrodynamical modeling of these systems. 

6. SUMMARY 

Our discovery of an abundant galaxy population at z ^ 
1.7 with extremely high emission line equivalent widths 
implies that many high-redshift, low-mass galaxies form 
many of their stars in extreme starbursts. We propose 
that we have observed an important formation mode for 
dwarf galaxies: a small number of strong starbursts that 
occur at early epochs {z > 1) each form ^ 10* Mq in 
stars in a very short time span (~ 30 Myr) to build up 
the bulk of the stellar components of present-day dwarf 
galaxies. This is in quantitative agreement with 'archaeo- 
logical' studies of present-day dwarf galaxies, which have 
shown that their star formation histories are burst-like 
and that the ages of their stella r populations sugg est for- 
mation redshifts z > 1 fe.g.. IWeisz et al1l2011t l. Un- 
der the reasonable assumption based on ACDM predic- 
tions for galaxy grwoth that the observed galaxies grow 
in mass by less than an order of magnitude up to the 
present day, our observations provide direct evidence for 
such an early formation epoch and, in particular, that 
short-lived bursts contribute much or even the majority 
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of star formation in dwarf galaxies. 

Facilities: HST(ACS,WFC3), 

Spitzer(IRAC) 



VLT(VIMOS), 
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